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Abstract 



This brief note is devoted to a study of genuine non-perturbative corrections to the Lan- 
dau gauge ghost-gluon vertex in terms of the non-vanishing dimension-two gluon condensate. 
We prove these corrections to give account of current SU(2) lattice data for the vertex with 
I different kinematical configurations in tire domain of intermediate momenta, roughly above 

2-3 GeV. We pay special attention to the kinematical limit which the bare vertex takes for its 
I tree-level expression at any perturbative order, according to the well-known Taylor theorem. 

00 ■ Based on our OPE analysis, we also present a simple model for the vertex, in acceptable 

CO I agreement with the lattice data also in the IR domain. 

o\ ' 
o 



H ■ 1 Introduction 



The infrared properties of the Landau-gauge QCD Green functions have motivated many studies 
in the last few years, mainly involving both lattice (see for instance Refs. [1-8]) and continuum 
approaches (see for instance Refs. [9-17] using Dyson-Schwinger equations (DSE), [18-21] using 
the refined Gribov-Zwanziger formalism or [22, 23] using the Curci-Ferrari model as an effective 
description). Most of the DSE analysis take advantage of approximating the ghost-gluon vertex 
by a constant when truncating the infinite tower of the relevant equations. To be more precise, 
only the behaviour of the involved transverse form factor needs to be approximated by a constant 
for the purpose of truncating the ghost propagator DSE (see also the recent paper [24]). However, 
this ghost propagator DSE (GPDSE) appears to be a cornerstone for the analysis of infrared DSE 
solutions for the gluon correlations functions: the two classes of solutions, namely "decoupling" 
and "scaling", depend on the regular or singular behaviour of the ghost dressing function at 
zero-momentum and on how this constrains the gluon behaviour via the GPDSE [11]. Most of 
the ammo for the approximation of the ghost-gluon vertex and GPDSE truncation is mainly 
provided by the Taylor theorem which is widely known as a non-renormalization one [25]. This 
claims that, in the particular kinematical configurations defined by a vanishing incoming ghost- 
momentum, no non-zero radiative correction survives for the Landau-gauge ghost-gluon vertex, 
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which takes thus its tree-level expression at any perturbative order [25]. This statement entails 
two important consequences: that (i) the bare ghost-gluon vertex is UV-finite for any kinematical 
configuration, and that (ii) in the specific MOM scheme where the renormalization point is taken 
with a vanishing incoming ghost momentum (which we named the Taylor scheme -T-scheme- 
in [26]) the renormalization constant for the ghost-gluon vertex^ is exactly 1. Then, one can use 
those two conclusions to write (schematically) the renormalized GPDSE as 

= 1 + arif^') J d\K{k,q)H\-'%q,k)F{q^) , (1) 

where ^ is the subtraction momentum for all the renormalized quantities, F is the ghost- 
propagator dressing function, K(k, q) is the part of the kernel which does not involve the ghost- 
gluon vertex, ax is the running coupling in the T-scheme^, and H^^'^'^ is the transverse form 
factor of the finite bare ghost-gluon vertex, defined by 

Tf''{-q,k;q-k) = -gor^^ (^q^H,{q, k) + {q - k)^H2iq, k)^ . (2) 

The kernel K corresponds to the gluon propagator, its transverse character in the Landau gauge 
allows to project out the i/2-form factor from expression (2). Therefore, we have chosen to call 
the surviving piece, i.e. Hi, the transverse form factor. Extending this nomenclature, we call H2 
the longitudinal vertex form factor. The kernel K can be either built by invoking lattice data 
(as done, for instance, in ref. [12]) or determined by solving the coupled system of Eq. (1) and 
the gluon propagator DSE. In both cases, a precise determination of the ghost dressing function 
from Eq. (1), to be confronted for instance with lattice data, requires a correspondingly precise 
knowledge of the transverse form factor Hi. Thus, the ghost-gluon vertex plays an important 
quantitative and qualitative role for the infrared DSE analysis of Yang- Mills QCD Green func- 
tions. In fact, with the usual bare ghost-gluon vertex approximation, namely that the bare 
transverse form factor is 1 over all the momentum range, it is readily appreciated from e.g. the 
results of [12, 30] that when the lattice data for the gluon propagator is fed into the GPDSE, 
there is a qualitative difference with the ghost propagator lattice data, unless one boosts by hand 
the input for the coupling constant at the renormalization point to an artificially large value. 
Another way-out [12] is to impose the bare transverse form factor to take a constant, clearly 
different from 1, to borrow the boost of the coupling input needed to describe lattice data'^. This 
is strongly hinting towards a missing piece in the ghost-gluon vertex, which, if stronger than a 
tree level one, will lift the output of the GPDSE up to the ghost lattice data (see refs. [31, 32]). 
Of course, before this can be precisely tested in practice, we need to cook up a reliable input for 
the ghost-gluon vertex. How to do this will be investigated in this work. 

On the other hand, the study of the running of the QCD Green functions triggered some lat- 
tice works [33-44], many of them aiming to determining the Aqcd parameter, which concluded 

^It can also be straightforwardly concluded that the ghost-gluon vertex renormalization constant, Zi, is exactly 
1 in the MS scheme. 

^The coupling in T-scheme has been profusely computed on the lattice and confronted to perturbative pre- 
dictions to estimate A^jg- [26-28]. It has been also used to define an effective charge that could be applied for 
phenomenological purposes [29]. 

^As it is shown in refs. [31, 32], taking a constant larger than 1 for the bare Hi helps to account better, but 
not enough, for the ghost propagator lattice data with the results of the integration of the GPDSE. Some non- 
perturbative structure for the ghost-gluon vertex is thus needed, and it is also preliminary shown there that the 
sort of inputs we can obtain with the OPE analysis in this work might very well lead to a good description of 
lattice data. 
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that the impact of the non-perturbative corrections coming from the dimension-two gluon con- 
densate {A'^) was not neghgible [37-41]. In particular, in ref. [26] we computed the running 
couphng renormahzed in the Taylor scheme from the lattice in pure Yang-Mills gauge theory 
and proved the necessity of including non-perturbative power corrections to describe properly 
the running over a momenta window from around 3 to 5 GeV. In the T-scheme, the configura- 
tion of momenta for the ghost-gluon vertex is the one "protected" by the non-renormalization 
theorem and, at least in perturbation theory to any order, no correction to the bare ghost-gluon 
vertex is expected beyond the tree-level expression. However, non-perturbative OPE corrections 
should be expected for any other kinematical configuration and, furthermore, it is interesting 
to examine whether a genuine non-perturbative correction might survive in the particular con- 
figuration where Taylor's argument works. More interesting will be to investigate whether, for 
the same price, these non-perturbative corrections related to {A^) might suggest, apart from 
describing the current lattice data inside the appropriate momenta window, a reliable input for 
the ghost-gluon vertex we mentioned above. This is the main purpose of this note. 

2 The OPE for the ghost-gluon vertex 

The OPE procedure for the Landau-gauge ghost-gluon vertex, 

V^'%-q,k;q- k) = V'^}'''^' {-q,k-q - k) G%,{q - k) F-'(g) F'^' {k) 

= I d^y d\ e*(«-'=)-^ e'^'-y ( T (c^(y)4(x)c'^(0)) ) (3) 

is quite similar to the one described in ref. [26]. Here, the OPE expansion shall read 
V^'^{-q,k-q-k) = ido)f'{q,k) 

\abcfi'i/' 



+ {d^rp^/ (Q^k) ^:A^,(OK'(0):) + ■■■ (4) 
where do accounts for the purely perturbative contribution to the vertex, while 

= 21^ + 2/W + 2/Pl + 4/[3] + 2/W + (5) 



where^ 



q-k N 




/W = 

- 2 q'^{q-k)'^ '^troe.Atl q,'^,q 1^) ■, {^) 
the tree-level ghost-gluon vertex being 

V^ttA-q, k-q-k) = -i p^2(g_fc)2 f"'" 9^,^^ - k) , (7) 



Nc is the number of colours. 



while 



and 




The other contributions, /['^'^'^l, can be immediately written using the previous results: 



and 




In all cases, the blue bubble refers to a contraction of the colour and Lorentz indices with 
^^abdfiu- We have also introduced the notation gj[y{t) = g^^u — for the transverse projector. 

Then, one obtains 

w-^' = {^y{q^k)-,SF{k)-,SF{q)^SG{k-<i)\^ V,tl^{-q,k;q - k) (12) 

with 

Nc 

saiq) = SF{q) = , 

I ,N Nc {q-k)-q [k - q) ■ k k ■ q\ 
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sp (resp. sg) comes from the OPE corrections to the external ghost (resp. gluon) propagator, 
i.e. from the non-proper diagrams in Eqs. (10,11) and sy from the proper vertex correction in 
Eqs. (6,9). The OPE corrections to the proper vertex can be obtained from there, remembering 
that the bare ghost-gluon vertex reads as: 

r-f"^^-q,k;q-k) = -g^P^' [q^H^{q,k) + {q - k) ^H^iq^k)^ . (14) 

Then, in Landau gauge, only the form factor Hi survives in the Green function to give: 

V^''-{-q,k-q-k) = -igor'%,gi,^{q-k)Hi{q,k)G{{q-kf)F{q^)F{k^), (15) 

where G and F are the gluon and ghost dressing functions for which we previously computed 
the non-perturbative OPE corrections. Thus, one would have: 

Hiiq,k) = Hr\q,k) (^l + sviq,k) 4(j^f ! + Oig\q-\k-\q-'k~^) ) . (16) 

All the non-proper corrections have been removed in the usual way. 

3 A model for the Euclidean ghost-gluon vertex 

We will now devote this section to model the transverse form factor of the Euclidean ghost-gluon 
vertex. Hi, on the ground provided by Eq. (16), and compare the result with some recent lattice 
data for this form factor computed in different kinematical configurations. 

3.1 The model 

In Euclidean metrics the bare ghost-gluon vertex can be written very similarly as: 

n'^reJ-Q, k-q-k) = igor'' (^q^^Hiiq, k) + {q - k),,H2{q, k)^ , (17) 

where the form factor Hi plays a crucial role when solving the ghost-propagator Dyson-Schwinger 
equation (GPDSE) as discussed above. The OPE non-perturbative corrections to the form factor 
Hi were obtained in Eq. (16), but that result is in principle only reliable for large enough k,q 
and q — k since the SVZ factorization on which it relies might not be valid for low momenta. 
Nevertheless, in the following we will propose a very simple conjecture to extend Eq. (16) to 
any momenta, which will provide us with a calculational model for the ghost-gluon vertex to 
continue research with. In particular, the perturbative part of the ghost-gluon vertex is usually 
approximated by a constant behaviour^; if we then apply a finite renormalization prescription 
such that 

Zi{^^)Hiiq,k)\^, = 1 , (18) 

where the renormalization momentum, fi'^, for a given kinematical configuration (for instance, 
q — k = and q^ = k"^ = /i^) is chosen to be large enough, on the basis of Eq. (16), we can 

^At least for large momenta, this seems to be the case in lattice simulations [45, 46] for several kinematical 
configurations, and so is confirmed also by the perturbative calculations in Refs [47,48]. 
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conjecture that 




1 + 



Ncg^ {A^) 



X 



^k'^cp- cos Q 



+ 2 



— ■s/Wcp' cos Q 



+ k"^ — 2\J q^k'^ cos 0) + mjj 



IR 



+ 2 



P(g2 + A;2 - 2v/^cos 0) + mfj 



/c^ — -s/W^ COS 6* 



IR 



) 



(19) 



gives a reasonable description of the ghost-gluon form factor Hi all over the range of its momenta 
q and k, where 9 stands for the angle between them. The purpose of Eq. (19) is to keep the 
main features of the momentum behaviour of the ghost-gluon form factor provided by the OPE 
analysis and to give a rough description of the deep IR only with the introduction of some IR 
mass scale, miR, which is mainly aimed to avoid the spurious singularities resulting from the 
OPE expansion in momentum inverse powers. 

It should be noted that the OPE technology implies first a hierarchic expansion in terms of 
local operators and next a factorization assumption that allows for all the large-distance infor- 
mation to be encoded in the condensates of those local operators in the QCD vacuum. If the 
momentum to be described decreases, more condensates of higher and higher dimensions should 
be needed, but also the Wilson coefficient containing the short-distance information should be 
computed at a higher loop level. To try to sum up all this is, if not impossible, a daunting 
task. This is why the OPE is applied to provide only an asymptotic expansion which is reliable 
to account for some QCD phenomenology, but nevertheless it can also successfully reproduce 
quenched and unquenched lattice results for many Green functions in some intermediate mo- 
mentum region [26,39-41]. Then, it is important to stress here that the model vertex (19) 
will never be obtainable from the OPE in QCD, as it is just impossible to describe with the 
OPE whatever quantity over a full momentum range. Very recently, a Landau gauge version 
of the Curci-Ferrari model has been shown to give, at the first perturbative order, a reliable 
description of quenched QCD lattice results [22, 23]. There is no reason to think that the same 
approach could also account for the operator expansion and factorization sum rules. As a matter 
of fact, one would obtain that applying the OPE technology with a massive gluon propagator, 
for instance with a frozen mass M, implies to replace 



in Eqs. (13,16) and this does not totally cure the above mentioned spurious singularities coming 
from the expansion in momentum inverse powers. In fact, assuming that M is the typical mass 
scale of a theory, the OPE machinery only intends to describe physics for momenta (much) larger 
than M. However, lattice results, as can be seen in the next subsection, and the Taylor theorem 
in its particular limit prove that no singularity appear in the IR. This is why we introduce in 
Eq. (19) the mass scale, mjR, that can be thought as some sort of IR regulator that should be 
related to some kind of effective gluon mass scale. 

The proposed expression is thus supposed to describe the vertex for the whole range of 
possible momenta, as needed to be plugged into the GPDSE in future work, while any OPE 
analysis, from its very definition, can only be valid for relatively large momenta. The proposal 
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(20) 



{q - kY 



{q - kf + M2 
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(19) has however the right asymptotic momentum behaviour, as compared to the OPE result 
(16). Its vahdity for lower momenta will shortly be checked by a confrontation with the available 
lattice data in the next subsection. Said otherwise, Eq. (19) can be also seen as a parametrization 
to describe lattice data (to be properly plugged into the GPDSE) that, behaving asymptotically 
as the OPE predicts, benefits of the fact that free parameters as {A"^) can be borrowed from 
(or compared with) different and independent analyses. It is anyhow very important to stress 
that it is the OPE analysis leading to (16) providing us with the kinematical structure of 
(19), at least in the intermediate region, and that the infrared completion has nothing to do 
with this. Thus, the success in the confrontation with the available lattice data, with different 
kinematical combinations for the form factor, will also confirm that the OPE analysis works in 
the appropriate region and it also gives an independent estimate of {A'^). 

Furthermore, to give further credibility to the vertex (19), we also plan a future study of the 
ghost-gluon vertex from the analytical viewpoint. For example, it would be instructive to find 
out if the refined Gribov-Zwanziger (RGZ) computational scheme could give access to a vertex of 
the type modeled here, since this formalism incorporates effects of non-perturbative dimension 2 
condensates in a controllable, renormalizable setting, applicable to the whole range of momenta. 
It is worth to remark that, as an enlightening indication, the gluon propagator in the RGZ 
scheme can be asymptotically expanded to agree with its OPE prediction stemming from the 
gluon condensate, {A'^); however, the deep infrared structure of the propagator clearly escapes 
from the OPE analysis and depends, for instance, on additional dimension two-condensates and 
on the Gribov-Zwanziger parameter defined through the so-called Gribov horizon condition [18- 
21]. An alternative could be the above-mentioned Landau gauge version of the Curci- Ferrari 
model, which was already tested at the propagator level in [22, 23] but that can be also applied 
to give a prediction for the ghost-gluon vertex. 

Of course, /.i^ being large enough, Zi can be approximated by some constant value for 
any /u^ because the logarithmic behaviour of the ghost-gluon in perturbation theory has been 
proven to be very smooth [47,48]. Thus Eq. (19) provides us with a very economical model 
because one needs nothing but some infrared mass parameter which, being related to the gluon 
mass, is expected to be 1 GeV, to parametrize the deep infrared behaviour of the ghost- 
gluon transverse form factor. It should be noted that, although the ghost-gluon vertex form 
factors do not diverge, the gluon condensate does: it is the product of the Wilson coefficient 
and of the condensate which is expected to remain finite thanks to a delicate compensation of 
singularities, in the sum-rules approach. Then, both the Wilson coefficient and the condensate 
should be renormalized, by fixing some particular prescription, and both shall depend on a 
chosen renormalization momentum. This is explained in detail in the work of ref. [49] where the 
Wilson coefficient for the quark propagator is included at order 0{a^). Since we carry out the 
computation of the Wilson coefficient in Eqs. (6-11) only at tree- level and neglect its logarithmic 
dependence, we do not specify any renormalization momentum for the gluon condensate. 

3.2 Comparing with available lattice data 

Some SU(2) and SU(3) lattice results for the ghost-gluon vertex are available in the literature 
(see for instance Refs. [45,46]). In particular, Cucchieri et al. have computed the tree-level- 
tensor form factor in SU(2) for three kinematics {p = q — k is the gluon momentum and ip the 
angle between the gluon and ghost momenta) : (1) = and if = tt/2, (2) p^ = and (3) 
If = 7t/3 with p^ = q^ = k'^ ([45]). As can be seen in the left pane of Fig. 1, when we plug 
g^{A^) = 6 GeV^, miR = 1.4 GeV and ^ = 1.04 into Eq. (19), the prediction for the form 
factor Hi appears to agree pretty well with the SU(2) lattice data for the three kinematical 
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configurations above mentioned. It should be also noted that the OPE prediction without the 
infrared completion through the introduction of miR given by Eqs. (13,16), plotted with dotted 
lines, also accounts very well for lattice data in the intermediate momenta region, above roughly 
2.5 GeV. Furthermore, the value for the gluon condensate lies in the same ballpark as the 
estimates^ obtained from the SU(3) analysis of the Taylor coupling in ref. [26]. This strongly 
supports that the OPE analysis indeed captures the kinematical structure for the form factor 
Hi. On the other hand, the infrared mass scale, as supposed, is of the order of 1 GeV and 
the (very close to 1) value for Z^^ accounts reasonably for perturbative value of Hi that we 
approximated by a constant. Thus, after paying the economical price of incorporating only 
these two more parameters, both taking also reasonable values, we are left with a reliable closed 
formula for the form factor at any gluon or ghost momenta. This is a very useful ingredient for 
the numerical integration of the GPDSE. 
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Figure 1: (left) The predictions of the model for Hi (solid lines) in eq. (19) and the OPE result without 
infrared completion (dotted lines) in eq. (16) confronted to the SU(2) lattice data borrowed from ref. [45] 
for the three kinematical configurations described in the main text: (fc — q)"^ = for ip — n/3 (violet) 
and (fi ~ tt/2 (red) and {k — q)^ ~ (green), (right) The results of the model for Hi in the SU(3) case 
and for the following kinematical configurations: = q^ for = 7r/3 (violet) and 9 = tt/2 (red) and 
— (green). 

The SU(3) results published by the authors of ref. [46] appear to be very noisy and cannot be 
invoked to properly discriminate whether a constant behaviour close to 1 or Eq. (19) accounts 
better for them. However, we can now take the mass parameter, mm, from the previous SU(2) 
analysis and the well-known SU(3) value for g'^ (^^), assume Z^^ = 1 and use Eq. (19) to predict 
the ghost-gluon transverse form factor, Hi. This is shown in the right plot of in Fig. 1 for three 
different kinematical configurations and, as can be seen, the deviations from 1 appear to be 
very small in all the cases and compatible with the results shown in Fig. 4 of ref. [46] for the 
vanishing gluon momentum case. 

4 The Taylor kinematics and the asymmetric gluon-ghost vertex 

An especially interesting kinematical configuration in the ghost-gluon vertex is the one where 
the incoming ghost momentum goes to zero {k — )> 0). Let us pay some additional attention 

®In ref. [26], g'^{A'^) is evaluated through an OPE formula with a tree-level Wilson coefficient to be 5.1 GeV^ 
at a renormalization point of 10 GeV. 
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to it. This particular kinematical limit is the one where Taylor's non-renormalization theorem 
provides us with a result for the bare vertex, which is at least exact up to all perturbative^ 
orders: 

n^roA-Q,0;q) = -gr'^q,.- (21) 

At any order of perturbation theory, this implies that Hi(q,0) + H2{q,0) = 1. According to 
Eq. (19), one would have 

Hi{q,0) = (l + Nc fJf\. ) • (22) 

V 4(iV^ - 1) g4 + J 

An interesting question to investigate is whether such a genuine non-perturbative correction still 
survives for the full ghost-gluon vertex in the Taylor limit, and not only for the transverse form 
factor Hi. Said otherwise, does eq.(21) maintains its validity upon adding the OPE corrections 
related to {A^) to it, or not? 

In order to properly address this question, one does not to have to study the ghost-gluon 
vertex for asymptotically large legs' momenta, as done in sec. 2, but to study instead the 
asymmetric ghost-gluon vertex 

V^''^{-q,e;q-e) = Vji,'''^' {-q,e;q - e) G%,{q - e) F'^'^' (q) F'^'^' (e) 

= J (fx e^^-?-")-^ ^ T (?(e)A|^(x)c'^(0)) ) (23) 

where, when taking the limit 0, c"(e) stands for an incoming vanishing ghost field. A certain 
care is needed since we have to deal here with a "so/it" ghost leg, similarly to the situation in 
ref. [40] when analyzing the three-gluon vertex. According to this, the key point for the OPE 
analysis of the asymmetric vertex is to apply the OPE expansion only to the fields carrying 
hard momenta and take then advantage of the operator expressions to compute the appropriate 
matrix elements in presence of a ghost field with "soft" momentum. Thus, we can apply to the 
asymmetric ghost-gluon vertex a procedure similar to the one outlined in sec. 2 and obtain 

V^'^i-q^e^q-e) = | d4xe*(''-^)- (^^(e) T (4(x)c«(0)) ) 

= K)?c'(9) :c^'(0):) 

+ (d,)ffX'iq) (^^(^) :c^'(0)<,(0)A^:(0):) + ••• (24) 

Then, if one follows standard SVZ techniques to obtain the perturbative expansion of OPE 
Wilson coefficients, we will take the perturbative vacuum for the matrix elements of Eq. (24) 
and obtain: 

V^t\A-Q, e-q-e) = {di)% {q, e) F-'(e) , (25) 

where di can be straightforwardly identified with the perturbative expansion of the ghost-gluon 
Green function, V, with an amputated soft-ghost leg. Beyond this purely perturbative first 

^The same argument of Taylor's perturbative proof still works if one consider the Landau-gauge ghost-gluon 
vertex DSE: a vanishing ghost momentum entering in the vertex implies the contraction of the gluon-momentum 
transversal projector with the gluon momentum itself for any dressed diagram. Thus, Taylor's theorem is still in 
order within the non-perturbative DSE framework. Some attention have been recently paid to the ghost-gluon 
vertex DSE [50]. 
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contribution, we will take the matrix element of Eq. (24) with two external soft gluons in the 
perturbative vacuum: 



id3)t'X' ^) ( ^mH:, m";,, (0) : c^' mi:, (OX, (O) : 



' flC'l 

\abii'v' 



b" . 



(26) 



On the other hand, when considering a non-perturbative vacuum for the matrix elements of the 
local operators in the OPE expansion of Eq. (24) and for a vanishing incoming ghost momentum 
e^, we will apply the factorization approximation (shown to work in ref. [40]) that follows: 



l?{e) :c^ (0)^^1,(0)4,(0): 



^0 

F^"'(e) 



?(.)c^'(0) 0) (:<,(0)4',(0) 



4(iV2 



(27) 



c 



Thus, Eq. (24) can be written as follows: 
V^''%-q,e;q-e) 



c 



I] 



F^'=(e) + 



(28) 



w 



abc 



where, now from Eq. (26), 

^ j^4^gig.(x-j;) hc^e)A''^,,{Q)Ai:,{id) T {A^^{x)c%y)) 



xa'b' 



(29) 



However, in the Landau gauge, the transverse incoming gluon kills the longitudinal contri- 
bution for the ghost-gluon vertex and only the transverse form factor. Hi, survives. One needs 
then to compute the proper (IPI) ghost-gluon vertex, which can be formally writen as 



K''{-q,k;q-k) 



Gfl,{q-k)F''^'{q)F^^'{k) ' 



(30) 



It is clear that, because of Eq. (25), one can factor the vanishing ghost propagator out in the 



r.h.s. of Eq. (28) and this implies that in Eq. (29) directly provides with the correction for 

the ghost-gluon vertex with the vanishing "soft" ghost leg explicitly amputated. Thus, we can 
formally define 



^fia) = vf'iQ) F'^'^iQ) G%{q) , 
such that the proper ghost-gluon vertex can be obtained as 



rf (-g,e;g-e) 



^abc 
pcrt,/^ 



{-q,£;q-£) + v\ 



abc 



4(iV2 - 1) 



+ 



where 



~abc 



2/[i] 



(31) 



(32) 



(33) 
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the factor 2 coming from the combinatorics and jt^l being the following diagram (no external 
legs) 



q-e 



/[I] 




Then, according to Eq. (2), one finally obtains: 

Tf''{-q,e;q-e) = -gf^'^ (q^ H^{q,e) + {q - e)^ H^iq^e)^ + 



where 



rPert. ^ ^ ^2 (Q ' ^) ' Q i^^) 



Thus, after taking the limit e — )• 0, one will have: 



(35) 



H,iq,e) = Hf {q,e) + Nag _ 
H.iq,s) ^ Hr\q,e) - No g^ • (36) 



rt!"'{-q,0;q) = -gP'' q^(^Hr\q,Q) + Hr\q,Q)^ = -gt'^ q^. . (37) 

This confirms that no non-perturbative OPE correction survives in the proper ghost-gluon vertex 
in the Taylor limit, although both form factors Hi and H2 separately undergo such a kind of 
correction. It is also straightforward to see that Hi{q,0) obtained from Eq. (36) would be 
the same as the one given in Eq. (22), obtained for the general kinematics, when q'^ ^ '^ir- 
Furthermore, provided we introduce in Eq. (36) for Hi the same infrared mass scale, mm, we 
used to obtain Eq. (22) we would again be left with the same result. These observations underline 
a rewarding agreement of our resuts although they stem from slightly different hypotheses. 



5 Conclusions 

The momentum behaviour of the ghost-gluon vertex has been studied in the framework of the 
operator product expansion (OPE). This approach has already been proved in the literature to 
be very fruitful in describing the running of two-point gluon and quark Green functions and of the 
strong coupling computed in several MOM renormalization schemes. In particular, the coupling 
derived from the ghost-gluon vertex renormalized in T-scheme, which can be directly computed 
from nothing else but the ghost and gluon propagators, led to a very accurate determination of 
Aqcd when its running with momenta obtained from the lattice estimate has been confronted 
with the OPE prediction, although only after accounting properly for a dimension-two gluon 
condensate, {A"^). In this note, we followed the same approach to study the ghost-gluon vertex 
and payed special attention to the transverse form factor for the ghost-gluon vertex, precisely 
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the one that plays a crucial role for the truncation and resolution of the GPDSE and that, 
in this context, is usually approximated by a constant. In this way we obtained a genuine 
non-perturbative OPE correction to the perturbative part for the transverse form factor. A 
very simple conjecture was then made to extend the OPE description beyond the momentum 
range where the SVZ factorization is supposed to work and we were left with a simple model 
describing the momentum behaviour of the ghost-gluon form factor. This model was shown to 
agree pretty well with a lattice SU(2) computation for several kinematical configurations of the 
ghost-gluon vertex, when the value for the gluon condensate is in the same ballpark as the one 
that is estimated from the running of the T-scheme coupling in SU(3) lattice gauge theory. This 
successful comparison with the available lattice data provides us with strong indications that (i) 
the OPE framework helps to account for the kinematical structure of the ghost-gluon vertex at 
intermediate momentum domain and (ii) that it helps to cook up a reliable model, providing a 
good parameterization also for the IR domain, to be plugged into the GPDSE to reproduce the 
ghost propagator lattice data. These are the main goals of the paper. 

Finally, we also proved that, in the particular kinematical limit associated with the Taylor 
theorem and hence by the T-scheme (a vanishing incoming ghost momentum), the corrections 
to the longitudinal and transverse form factors cancel against each other and that the tree-level 
result for the ghost-gluon vertex is thus recovered. 

In future work, we foresee to put to use the vertex (19) in combination with the GPDSE, this 
in order to investigate at the quantitative level how this vertex allows to bring the analytical 
predictions of the GPDSE in closer agreement with the lattice output for the ghost propagator, 
given the input of the kernel via the lattice gluon propagator. 
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